265 research outputs found
Quantifying admissible undersampling for sparsity-exploiting iterative image reconstruction in X-ray CT
Iterative image reconstruction (IIR) with sparsity-exploiting methods, such
as total variation (TV) minimization, investigated in compressive sensing (CS)
claim potentially large reductions in sampling requirements. Quantifying this
claim for computed tomography (CT) is non-trivial, because both full sampling
in the discrete-to-discrete imaging model and the reduction in sampling
admitted by sparsity-exploiting methods are ill-defined. The present article
proposes definitions of full sampling by introducing four sufficient-sampling
conditions (SSCs). The SSCs are based on the condition number of the system
matrix of a linear imaging model and address invertibility and stability. In
the example application of breast CT, the SSCs are used as reference points of
full sampling for quantifying the undersampling admitted by reconstruction
through TV-minimization. In numerical simulations, factors affecting admissible
undersampling are studied. Differences between few-view and few-detector bin
reconstruction as well as a relation between object sparsity and admitted
undersampling are quantified.Comment: Revised version that was submitted to IEEE Transactions on Medical
Imaging on 8/16/201
High resolution image reconstruction with constrained, total-variation minimization
This work is concerned with applying iterative image reconstruction, based on
constrained total-variation minimization, to low-intensity X-ray CT systems
that have a high sampling rate. Such systems pose a challenge for iterative
image reconstruction, because a very fine image grid is needed to realize the
resolution inherent in such scanners. These image arrays lead to
under-determined imaging models whose inversion is unstable and can result in
undesirable artifacts and noise patterns. There are many possibilities to
stabilize the imaging model, and this work proposes a method which may have an
advantage in terms of algorithm efficiency. The proposed method introduces
additional constraints in the optimization problem; these constraints set to
zero high spatial frequency components which are beyond the sensing capability
of the detector. The method is demonstrated with an actual CT data set and
compared with another method based on projection up-sampling.Comment: This manuscript appears in the proceedings of the 2010 IEEE medical
imaging conferenc
An algorithm for constrained one-step inversion of spectral CT data
We develop a primal-dual algorithm that allows for one-step inversion of
spectral CT transmission photon counts data to a basis map decomposition. The
algorithm allows for image constraints to be enforced on the basis maps during
the inversion. The derivation of the algorithm makes use of a local upper
bounding quadratic approximation to generate descent steps for non-convex
spectral CT data discrepancy terms, combined with a new convex-concave
optimization algorithm. Convergence of the algorithm is demonstrated on
simulated spectral CT data. Simulations with noise and anthropomorphic phantoms
show examples of how to employ the constrained one-step algorithm for spectral
CT data.Comment: Submitted to Physics in Medicine and Biolog
Estimating the Spectrum in Computed Tomography Via Kullback–Leibler Divergence Constrained Optimization
Purpose
We study the problem of spectrum estimation from transmission data of a known phantom. The goal is to reconstruct an x‐ray spectrum that can accurately model the x‐ray transmission curves and reflects a realistic shape of the typical energy spectra of the CT system. Methods
Spectrum estimation is posed as an optimization problem with x‐ray spectrum as unknown variables, and a Kullback–Leibler (KL)‐divergence constraint is employed to incorporate prior knowledge of the spectrum and enhance numerical stability of the estimation process. The formulated constrained optimization problem is convex and can be solved efficiently by use of the exponentiated‐gradient (EG) algorithm. We demonstrate the effectiveness of the proposed approach on the simulated and experimental data. The comparison to the expectation–maximization (EM) method is also discussed. Results
In simulations, the proposed algorithm is seen to yield x‐ray spectra that closely match the ground truth and represent the attenuation process of x‐ray photons in materials, both included and not included in the estimation process. In experiments, the calculated transmission curve is in good agreement with the measured transmission curve, and the estimated spectra exhibits physically realistic looking shapes. The results further show the comparable performance between the proposed optimization‐based approach and EM. Conclusions
Our formulation of a constrained optimization provides an interpretable and flexible framework for spectrum estimation. Moreover, a KL‐divergence constraint can include a prior spectrum and appears to capture important features of x‐ray spectrum, allowing accurate and robust estimation of x‐ray spectrum in CT imaging
Toward optimal X-ray flux utilization in breast CT
A realistic computer-simulation of a breast computed tomography (CT) system
and subject is constructed. The model is used to investigate the optimal number
of views for the scan given a fixed total X-ray fluence. The reconstruction
algorithm is based on accurate solution to a constrained, TV-minimization
problem, which has received much interest recently for sparse-view CT data.Comment: accepted to the 11th International Meeting on Fully Three-Dimensional
Image Reconstruction in Radiology and Nuclear Medicine 201
- …